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INTRODUCTION 


In  this  report  we  display  the  results  of  analysis 
of  microbarograph  records  from  two  presumed  atmospheric 
nuclear  explosions.  The  events  occurred  on  27  December 
1968  and  29  September  1969,  and  were  recorded  by  micm- 
barograph  arrays  in  Boulder,  Colorado  (BC) ,  College, 

Alaska  (CA)  ,  lluancayo ,  Peru  (IIP),  Tel  Aviv,  Israel  (TA) , 
and  at  LAMA,  the  Large  Aperture  Microbarograph  Array  in 
Montana.  We  also  show  LAMA  data  and  analysis  results 
for  eight  other  presumed  explosions.  For  array  coordinates 
and  locations,  see  Appendix  I. 

In  addition  to  demonstrating  the  effectiveness  of 
our  digital  detector,  the  pupose  of  this  report  is  to 
make  the  calibrated  data  and  our  spectral  calculations 
available  to  other  workers  in  this  field. 

The  LAMA  data  were  recorded  digitally,  and 
calibrations  were  available  for  the  two  dates  of  the 
events  studied  here.  However,  the  other  data  recorded 
on  analog  magnetic  tape,  and  severe  difficulty  was 
experienced  in  locating  and  processing  the  signals  from 
these  events.  First,  a  time  control  channel  is  included 
among  the  analog  tape  channels,  but  the  hour  marks  were 
not  identified,  nor  were  the  days.  Each  tape  contains 
many  days  of  data,  so  that  locating  the  desired  event 
requires  that  the  hours  be  counted  from  some  identifi¬ 
able  point  on  the  tape.  After  long  and  frustrating 
attempts  to  establish  a  time  base,  we  finally  had  to 
rely  on  the  time  notations  written  at  infrequent  and 
irregular  intervals  in  marker  pen  ink  on  the  back  of 
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the  analog  tape  by  the  station  operators. 

Counting  time  marks  from  such  benchmarks,  however, 
assumes  continuity  of  recording,  and  the  records  from 
the  stations  recording  on  analog  tape  are  very  frequently 
intermittent.  The  fact  that  often  the  interruptions  were 
not  noted  on  the  station-  logs  made  it  impossible  to 
locate  signals  on  the  tapes.  The  results  presented  here 
are  those  for  station  records  which  we  were  eventually 
able  to  locate  and  identify. 


DATA  I' K  LI*  A  RAT  I  ON 


A  great  amount  ot  *‘ffort  was  required  simply  to 
prepare  th-’*  data  fro*  outlying  stations  for  digital 
analysis.  We  will  go  into  some  detail  about  the 
procedures  we  used,  since  the  unexpected  delay  was  the 
major  reason  why  the  project  ended  before  meaningful 
coclusions  could  he  drawn  frc<:?  t.ic  data  analysis. 

A  number  of  steps  were  involved: 

1.  Locate  events  and  calibrations  on  AM  tapes.  The 
timing  channel  on  the  analog  tapes  does  not  distinguish 
the  time  marks;  every  hour  mark  looks  the  same  as  every 
other  hour  mark.  The  tapes  run  for  weeks  at  a  very  slow 
speed,  and  the  station  operators  occasionally  mark  the 
day  and  hour  on  the  record  with  a  felt-tip  pen,  often 

at  intervals  of  many  days  or  weeks.  Since  the  signal-to- 
noisc  ratio  of  the  events  being  sought  was  often  so  low 
that  the  signals  were  not  visible  on  the  single  channels, 
we  had  to  rely  on  the  timing  information.  However,  the 
stations  frequently  stopped  operation  and  started  again 
after  an  outage  of  hours  or  d»ys,  and  it  was  not  uncommon 
that  the  station  personnel  did  not  note  the  outage  in 
the  station  logs.  The  signals  reported  here  were 
finally  located  with  great  difficulty. 

2.  Convert  from  AM  analog  tape  to  TM  analog  tape. 

The  digitizing  equipment  at  the  Seismic  Data  Laboratory 
did  not  have  the  proper  characteristics  to  digitize  the 
AM  analog  ttpes,  so  wc  had  to  use  an  intermediate  step, 

in  which  the  AM  tape  was  dubbed  onto  an  FM  tape  at  30  ips. 
This  shxitcd  the  recorded  spectrum  up  from  0,02  cps 
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(low  end)  to  14-1  cps,  and  1  row  1.0  cps  (high  enJ)  to 
7200  cps,  which  the  digitizing  equipment  could  handle. 

Khen  they  could  he  located,  calibrations  were  also 
dubbed  onto  the  I'M  tapes. 

J.  A/0  conversion  of  FM  tapes.  After  applying  an 
anti-aliasing  low-pass  filter,  the  FM  tapes  were 
digit  l*.  cd  at  a  rate  which  corresponds  to  two  samples 
per  second  on  the  original  frequency  scale. 

4.  Convert  digital  tape  to  SOI.  library  tape  iorraat. 
This  essentially  involved  demultiplexing  the  data  channels. 

5.  Convert  SOI.  library  tapes  to  Subset  format.  This 
step  put  the  data  in  a  form  accessible  to  our  Fortran 
programs. 

Programs  used  in  preparation  and  analysis  of  the 
data  are  described  by  Kerr  (1971). 
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Host  beams  were  computed  using  program  ONH  (Kerr, 
1971)  on  the  inner  five  elements  of  the  LAMA  array  (for 
details  sec  Flinn,  1972).  This  process  finds  by  an 
iterative  process  the  array  beam  havirg  the  largest 
F-statistic: 

(N-l)  Vb/N 

VT  -  VN' 

where  N  is  the  number  of  channels,  Vg  is  the  variance 
of  the  beam  sum,  and  V^,  is  the  total  variance  of  all 
channels  within  the  time  window  under  consideration. 

Raw  data  of  the  station  recordings  are  shown  in 
Figures  1  through  9.  The  best  beam  records  are  shown 
in  Figures  9  through  16.  Power  spectra  computed  for  the 
first  34  minutes  of  signal  in  these  beam  sums  are  shown 
in  Figures  17  through  22.  Time-varying  spectrograms 
(Cohen,  1970)  were  computed  using  a  two-minute  window 
with  50  percent  overlap  (Figures  23  through  28).  The 
signal  spectral  estimates  were  also  computed  by 
following  the  best  beam  frequency  by  frequency  througn 
f-k  space.  The  power  at  each  frequency  in  the  estimate 
was  taken  to  be  the  maximum  power  in  the  wavenumber 
plane  at  that  frequency.  Four  spectral  plots  were 
generated  for  each  signal:  power,  F  statistic,  phase 
velocity,  and  azimuth  (Figures  29  through  35).  See 
Smart  and  Flinn  (1971)  for  discussion.  The  one  exception 
here  is  the  record  for  Tel  Aviv  on  29  September  1969, 
which  had  only  three  operating  channels  that  day. 


ANALYSIS  OF  OTHER  LAMA  RECORDINGS 


Microbarograms  of  eight  other  presumed  atmospheric 
explosions  were  recorded  at  LAMA,  and  were  analyzed  by 
the  'best  beam'  processor  described  above.  The  raw  data 
are  shown  in  Figures  36  through  41,  and  the  spectra  in 
Figures  42  through  49. 
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*  Information  from  National  Ocean  Survey, 
Preliminary  Determination  of  Epicenters  reports 


SEISMOGRAM  NO.  249S8 


Figure  8.  Event  of  29  September  1969  recorded  at  LAMA,  Montana, 


Figure  9.  Event  of  27  Decemb 
Beam  velocity  and  azimuth  ar 


e  11.  Event  of  29  September  1969  recorded  at  College,  Alaska, 
velocity  and  azimuth  are  330  meter/second  and  310°. 
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re  13.  Event  of  27  December  1968  recorded  at  Tel  Aviv,  Israel 
velocity  and  azimuth  are  310  meter/second  and  62°. 


Figure  20.  Event  of  27  December  1968  recorded  at  Tel  Aviv,  Israel 
Power  spectrum  computed  starting  at  12  10. 


Figure  21.  Event  of  2S  September  1969  recorded  at  Tel  Aviv,  Israeli 
Power  spectrum  computed  starting  at  13  20 <, 
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Figure  24,  Time-varying  spectrum  of  the  event  of  29  September  1969 
recorded  at  Huancayo,  Peru. 
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Figure  25.  Time-varying  spectrum  of  the  event  of  27  December  1968 
recorded  at  Tel  Aviv,  Israel. 
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Figure  26,  Time-varying  spectrum  of  the  event  of  29  September  1969 
recorded  at  Tel  Aviv,  Israel. 


35 


0  .008  .011  .016  .021  .021  .030  .035  .040  .045  .050  .055 

FREQUENCY  (cps) 


Figure  27.  Time-varying  spectrum  of  the  event  of  27  December  1968 
recorded  at  LAMA,  Montana. 
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ime-varying  spectrum  of  the  event  of  29  September  1969 
LAMA,  Montana. 
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Figure  29.  Best-beam  spectral  estimates  for  the  event  of  27  December 
1968  recorded  at  Boulder,  Colorado.  From  top  to  bottom:  power, 
ai rival  azimuth,  phase  velocity,  and  F  statistic.  Figures  34  through 
39  are  arranged  in  the  same  way. 
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Figure  31.  Best-beam  spectra  for  the  event  of  29  September  1969 
recorded  at  College,  Alaska. 
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Figure  32.  Best-beam  spectra  for  the  event  of  29  September  1969 
recorded  at  Huancayo,  Peru. 
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Figure  33.  Best-beam  spectra  for  the  event  of  27  December  1968 
recorded  at  Tel  Aviv,  Israel. 
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Figure  38.  LAMA  data  channels  for  16  May  1970. 


Figure  39.  LAMA  data  channels  for  23  May  1970. 


Figure  40 „  LAMA  data  channels  for  4  July  1970. 


Figure  42.  Best-beam  spectra  for  LAMA  data,  0400-0508,  4  August  1968 


Figure  44.  Best-beam  spectra  for  LAMA  data,  0200-0308,  9  September  1968 


SEISMOGRAM  1450)43 


Figure  46.  Best-beam  spectra  for  LAMA  data,  0145-0253,  23  May  1970 


Figure  48.  Best-beam  spectra  for  LAMA  data,  0145-0253,  4  July  1970 
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